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CHAPTER I 


Introduction 


Let R be a commutative Euclidean domain with identity 1 
and with map d, which assumes integral non-negative values 
for~ ali nonzero a, ineR,ssuch’ that for any a, bin Rk, b #¥ 0, 
there exist gq and r in R for which a = gqeb + r and either fr 
= 0 or d(r) < d({b). Then M(n,R), the set of nxn matrices 
over R, forms a left Euclidean domain provided R is proper 
(i-e., if R is not a field and d(aeb) = d(a)ed(b) for every 
a, b of R) [Sanov; 1967]. If R is proper, then a_ suitable 
map for M(n,R) is d(det(A)), where A belongs to M(n,R) and 
det (A) is the determinant of A. A more general case is 
discussed by Brungs [Brungs; 1973], but it is not considered 


here. 


In M(n,R), the following theorem on left division holds 


true [Sanov; 1967] : 
Theorem 1.1. Let A, B belong to M(n,R). If det(B) # 0, then 
there exist 0, R in M(n,R) such that 

(1) either A= Be Q, 

(2) or A=BeQ+Rand0O < 4 (det (R)) 

< d(det(B)). 

The Euclidean algorithm can therefore be applied in an 
obvious way to prove the existence of a left greatest common 


divisor D (which is defined in the chapter II) of two 
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matrices A, B in M(n,R). In chapter II, we show that the 
application of the Euclidean algorithm to construct a left 
greatest common divisor D is nontrivial. A more efficient 
procedure for this construction involves first the 


computation of the Smith normal form of both A, B. 


Theorem 1.2. (Smith normal form) [Newman; 1972,pp.26- 
27] Let A belong to M(n,R). Then there exist two unimodular 


matrices U, V in M(n,R) such that 


A=UeS eV, 
where 
S = diag (S)¢S,e--+++7S_) 
<4 < ne-tTe = 5 See = = ; 
and s.| Say ToS -in's/ n=-trrand S41 So Oy ie 
rank (A) = cr. 


The theorem states that A is equivalent to a diagonal matrix 
S. If in addition U and V can be determined so that V = U-1, 


then A is said to be similar to S. 


As we shall see, the construction of a matrix D given 
the Smith normal form of A, B is a simple procedure. 
Furthermore, in the chapter V, we show that the Bradley's 
algorithm for computing the Smith normal form of matrices is 
inexpensive (relative to the direct application of the 
EBuclidean algorithm to A, B). Therefore, the procedure we 
recommend for computing a left greatest common divisor of 
two matrices requires essentially the construction of the 


Smith normal form. 


The Smith normal form of a matrix A in addition yields 


on inspection a decomposition of the matrix into its prime 
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factors. It is perhaps surprising to note that unlike the 
case for the integers and polynomials the prime 
decomposition of matrices can not be applied directly for 


finding a left greatest common divisor. 


The decomposition of a matrix into its prime factors is 
interesting for its own sake, and in chapters III and IV we 
deviate somewhat from the main theme to explore this subject 
further. Perhaps the highlight in these two chapters is the 
introduction of a normal prime matrix which leads to some 


strong uniqueness results. 


In concluding this chapter, we remark that for the sake 
of simplicity, in many of the results to follow, R is 
restricted to be the ring of integers Z. The map d is then 
be the absolute value function "Jj i". In most cases, 
however, it should be clear that the results given can 


eaSily be generalized to arbitrary Fuclidean domains. 
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CHAPTER If 


The Euclidean Algorithm for Matrix 


Rings over the Integers 


2-1 : Euclidean Algorithm for Matrix Rings 


——— —S a = a SSS 


In Sanov's paper an algorithm for division in the ring 
M(n,Z) is given. His algorithm is briefly summarized in 
section 2.2. AS is true for integer and polynomial rings, 
once division is defined, the Euclidean algorithm can be 
used to find a greatest commom divisor of two elements in 


M(n,Z). First, however, we give some definitions ;: 


Definition 2.1. If A, B are in M(n,Z), B # O (the zero 


——— 


matrix) then B left divides A if A = B © C for some matrix C 


in M(n,2Z). 


Definition 2.2. Let A, B be in M(n,Z), one of which is 


nonsingular. A left greatest common divisor of A and B, 
denoted by lgcd(A,8), is a matrix D in M(n,Z) such that D 
left divides both A and 8B, and furthermore if D' is any 
other matrix in M(n,Z) which left divides both A and B then 


D' left divides D. (The nonsSingularity condition is relaxed 


in chapter V). 


Theorem 2.1. Let a left greatest common divisor D of A 
and B be nonsingular. If D' is another left greatest common 


divisor of A and B then 
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where U is a unimodular matrix in M(n,Z). 
Proof : 
By definition, there exists matrices X and Y such 


that 


Since D is nonsingular, det(Y e X) must be the 
identity element in Z, and X, Y are therefore both 
unimodular matrices in M(n,Z). 


QeE.D. 


Given an arbitrary matrix A and a nonsingular matrix B 
in M(n,Z), this algorithm finds their greatest common 
divisor. 

Step 1 (Division) : 
| Set 
Ro <=—/ A god 5 
bY <== 45 
By <== PRs 
“Step 2 (Termination) : 
If B = 0, then terminate with A as the answer; else 


go to step 1. 


To show the validity of the algorithm we first state 
without proof that if a matrix A left divides both B and C 


then A left divides B e X + C e Y for any matrices X, Y in 
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First, left divide A by B_ getting, according to 
Sanov's division algorithm, a quotient ot and a 
remainder Ry such that’ *aAes= 2B oe % & R) with 0 < 
1) = 0 then Ri = 0 
(see Theorem 1.1) and B left divides A so that 


Jdet(R,)1 < [det(B)1. If det(R 


lgcd (A,B) = B. If det (R,) # 0, we divide B by 


Ry getting a quotient Q and remainder R., such that B 


2 


=) RB) * 05 + BR, with 0s [det (R5y1 < {det (R,)1. If 
det (R,) = 0 the procedure again terminates; whereas, 
if det(R ) # 0 we repeat to obtain Ry = R, ® Q, + 


R, with 0 < |det(R,)I < [det(R,) I. 


Eventually the process must terminate with a zero 
remainder since the decreasing sequence of 
nonnegative numbers 

[det (B){ > Idet(R,) {| > Idet(R,)| > .....- 
can be repeated at most |det(B)]| times. We therefore 
obtain the equations ;: 
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We now show that Rus the last nonzero remainder, 
is a left greatest common divisor of A and B. Since 


and R, left divides R the 


k k-1" k ke 
next to the last equation in (1) implies that RL left 


R; left divides: R 


divides R This process may be continued to _ show 


k-2° 


that Ri left divides A and B. 


On the other hand, if some other matrix R' left 
divides A and B, then it follows from the _ second 
equation in (1) that R' left divides Ro. Continuing 
this argument step by step, we finally have that R!? 
left divides R. - Thus R, is a left greatest common 
divisor of A and B, so that lgcd(A,B) = Ri 

Oake Dis 


2.2 : Sanov's Rational-Arithmetic Division Algorithm 


A 


Short description of the Sanov's division algorithn 


for two matrices A, B in M({n,Z) with B nonsingular follows : 


Step 1 


Step 2 


Step 3 


(Triangulation) : 

Perform elementary column operations on A and B to 
get, vtespectively, lower triangular matrices A‘ and 
B'. Let X be the unimodular matrix such that A' = A ® 
9 Ge 

(Inversion) : 

Find the inverse matrices (B')~!, X71, 
(Multiplication) ;: 

Compute the lower triangular matrix Q = (q..) = 


(B')-1t e At, 


Step 4 (All components of Q are integers) : 
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Step 5 


Step 6 


If all components of Q are integers then return 
with ER <-— O07 (inGthis case B leftidividés!a) . 

(Q has non-integral diagonal elements) : 

If some elements on the diagonal of Q are non- 
integral, then construct a lower triangular matrix R! 


with components (c'55) as follows ; 


i, if q;; is an integer, 
ria 
Qn — tgs, 1 qj, iS’ not “an integer. 
By aay HLF ie 32 
Then set 


Re <-=2 BY ie TRS) -e) X-* 
and return. 
(Q has integral diagonal elements but non-integral 
off-diagonal elements) : 
Let the first non-integer element be Se ie A matrix 


R'* with components (r'! is then constructed as 


io 
follows : 
(1) : Diagonal elements : 


0; 26 23 eats, 
io otherwise. 


(2) : Upper off-diagonal elements ; 


0, otherwise. 


(3) : Lower off-diagonal elements :; 
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i : Along the sth off-diagonal : 


0% LEGGE 3, 
C's4s5 = Qj+siv LEW lane ys 
Vasa 7 “Viesit otherwise. 


ii : Otherwise : 
0, Thijs hb. < Ss 
Cig = 
Iige 
Then set 
R <-- Bt e Rt e X-1 


and return. 


Section 2.3 : An Integer-Arithmetic Division Algorithm 


Sanov's division algorithm involves arithmethic 
operations on rational numbers. Computationally such 
operations are undesirable, since in order to minimize the 
growth of intermediate results, common factors must be 
removed. This requires the computation of the greatest 
common divisor of numbers which could be large in magnitude. 
In this section, therefore, a new division algorithm which 
requires only integral arithmetic operations is given. At 
this stage of the research, however, no claims on the 
computational superiority of the new algorithm shall be 
made. 


Integer-Arithmetic Division Algorithm : 


Given an arbitrary matrix A and a nonsingular matrix B 


in M(n,Z), this algorithm finds (A mod B). 
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Step 1 


Step 2 


Step 3 


Step 4 


Step 5 


Step 6 


10 


(Computation of the adjoint matrix) : 
Find the adjoint matrix, Bt, of the matrix B. 
(Multiplication) ;: 
Set GPpCFetoreerne 
(Smith normal form) : 
Find the Smith normal form Ss with diagonal 
components (s;), of the matrix C such that 
S S20 er Ce V. 
{Split matrix) : 
The matrix S is split into the sum of two diagonal 
matrices Q*', R* with components (q'.), (c'.), 
respectively. The components ae rs are constructed 
as follows ;: 
(1) 2: If det(B) divides Sis set 
tt; = det(B), 
q'. = (Ss. f adet (B)) = te 
(2) : If det(B) can not divide S,, set 
cla op TES. f/ det(B)4 e det(B), 
q’'; = ts; / det (B)4. 
(R*' has all elements equal to det(B)) : 
Tf all elements of R' are equal to det(B), then 
return with R <-- 0. 
(R' has some element less than det(B)) : 
If some element of R* is not egual to det(B), then 
set 
R <-- A- Be (U0 © Qf € ¥) 


and return. 


of the Integer-Arithmetic Division Algorithm : 
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By Theorem 1.2, there exist two unimodular matrices 
UO, ¥ in M(n,2Z) such that 

BtieoR ='U eS © V, 
where S, with diagnoal components Sie is the Smith 
normal form of the matrix Bt ® A. But for alli, 
there are integers Q'.s ci. such that 
Ss; = q*,¢ det (By) + r'y 
with 

O76 Ew < det(B). 
Let Q*' and R' be diagonal matrices with ith diagonal 
elements q'. and Cc‘. respectively. There are two 
cases to consider ; 
Ci 2) EE EY, = det (2) for ail i= ty... s chen 

Bree (ATS Ulens: evry 


= Ue (det(B) e Qf + R') e V 


i) 


det(By) © U e (Q' + IT) e V 
By left multiplying the above equation by B, we 
get 
A= Be (0 6. (OMe Epi ey). 
Thus, B left divides A. 
(2), 02 fee det(B) for some i, then 
0 < {det(R")| < [det(B)]". 
On the other hand, 


Bt e A = Ue (det(B) © Q' + Rt) e V 


det(B) ee 0 ® Qt e V 
+ Ue R' eV (2) 
with 


0 < tdet(U ¢ R' © V)t < [det(B) I”. (3) 


ee 
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If 


5) 
it 


Bee, 2) CUCS OF © ¥)iz 


AS rae (ae Ore Vink Rp (4) 
and we now claim that the matrix R satisfies 
the condition 

0 < {det (R)] < Jdet(B)]. 

By left multiplying equation (4) by the matrix 
B+ and comparing with the equation (2), it 
follows that 

Bites Roa Uverk®e°Vi. 

Thus, det(R) # O, and moreover by inequality 
(3) 


[det (R) | 


{det (U e R' e V){ / Jdet (Bt) | 
< [det (B)["/7 |det (B) |"-! 
= jdet(B)|. 


Q.E.D. 


In concluding this chapter an example is given which 
illustrates the integer-arithmetic division algorithn. 
Example : 


Given two matrices 


gc 04 
A= | 1 
i 4, 
a wees) 
B= | | 
es | eee 


We first obtain 
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eee 044 
Bt = | 
t —4 Tis 
and 
G2 7014 
Bt © A = | | 
t3 44. 


By applying some algorithm for obtaining the Smith 
normal form, we get 
Bt e@ A 
Che) ab aie: Te 0 aie ie os 
=o tel i | | 
L3 -14 60 8440 14. 


That is, the Smith normal form S is 


ct 0,4 
Ses tl | 

CO: |. ae 

and 

r 2 -14 
U = | | 

Lif Sts, 

r 1 4 4 
v= | I 

t 0 14. 


Now split the matrix S according to the step 4 to 


obtain 
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The quotient Q is then given by 


cr One rar 


Note that 


1 = |det(R)| < Idet(B)! = 2. 
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CHAPTER III 


One-Side Decomposition 


for Matrix Rings 


3.1 : Definition of A Prime Matrix 


From Theorem 1.1 we know that the matrix ring M(n,R) 
over a Euclidean domain R is a left Euclidean domain. One of 
the most interesting properties of principle ideal domains 
(and therefore of Euclidean domains as well) is that such 
rings admit a theory of unique factorization {MacLane; 1967, 
pp- 154-155]. One method for decomposing a matrix into its 
prime factors is discussed by Sanov [Sanov; 1967]. In this 
chapter a different decomposition of a matrix in M(n,Z) into 
its prime factors is introduced. This decomposition permits 
us to exhibit certain uniqueness result which are not 


possible using Sanov's decomposition. 


With the demonstration of this fact in mind, we first 
give the following definition ; 
Definition 3.1. A matrix P in M(n,Z) is called a prime 
matrix (i.e., P is a prime element in M(n,Z)) if |det{(P)| > 
1 and P has no other left divisors besides unimodular 
matrices and matrices which are right equivalent to P 
(matrix B is said to be right equivalent to A if B= Ae 
VY for some unimodular matrix V). 


This definition is equivalent to saying that for any 
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decomposition, P = A e@ B, SteaPeaxeither WAt“or Briis 


unimodular. 


Theorem 3.1. The determinant of any prime matrix in 


M(n,Z) is a prime element in 2%. Conversely, if the 
determinant of some matrix in M(n,Z) is a prime element in 
Z, then it must be a prime matrix in M(n,Z). 
Proof : 
The second part of the theorem is obvious, and the 
proof of the first part can be found in [Sanov; 
1967). 


QED. 


From this theorem, we know that diag(1,...,p,1,...1), 


with p a prime in Z, iS a prime matrix immediately. 


A = Py @ P, Oe aa) OP (1) 
be a decomposition of a matrix A in M(n,Z) into the prime 
factors P,, 18 k < m. The existence of such a decomposition 
is proven by Sanov. Given the decomposition (1) and 
arbitrary unimodular matrices U ,..., U e LUcES clear that 
A = (P, © Ujt) © (U, © P, © UFZ) © .... 
Wy eee Ves mM ep Shad 


me m 


is another decomposition of A, Since, in addition, M(n,Z) is 
non-commutative, the question of uniqueness of the 
decomposition is a nontrivial one. To make this problem more 
tractable, Sanov restricts the prime matrices to be of 
lower-triangular type. This leads to unigueness results 


which are too confining. In the remainder of this chapter, 
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we discuss the decomposition of a matrix into normal prime 


factors. Some interesting results on uniqueness then follow. 


3.2 : Normal Prime Matrices 


— SSS SS SS SS eS ee 


Definition 3.2. A prime matrix P in M(n,Z) is a normal 


prime matrix if it is similar to diag{1, ..., 1,p), where p 


is a prime element in Z. 


Examples : 
1 3: Let 
ean 12 
ak | 


Because det(P) = 2 is a prime of Z, the matrix P is a 
prime matrix. Moreover, by taking 
r 1 34 


u= | | 


then 


Thus, P is also a normal prime matrix. 


212° Lee 


Because det(Q) = 3 is a prime of Z, the matrix Q is a 
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prime. On other hand, there does not exist any 
unimodular matrix U such that 
Ue Q e U-t = diag(1,3). 


Thus, the matrix Q is not normal prime. 


Definition 3.3. A decomposition of a nonsingular matrix 


A, which is not unimodular, in M(n,Z) into the form 


where Pie 1< k << mM, are commuting normal prime matrices, 
and U is a unimodular matrix is called a left normal 


Se SSS SS 


decomposition of A. 


3.3 : Existence of a Left Normal Decomposition 


From Theorem 1.2, we know that for any nonsingular 
matrix A in M(n,Z), A is equivalent to its Smith normal 
form diag(S,,SyreeeerS)e 1.e., there are two unimodular 
matrices U, Vin M(n,Z), Such that 3; 

A= Ue diag(S,+Sore+eeeS_) e Vv. (2) 
By standard arguments, we know that all components S,, 1 < k 
< n, can be decomposed into products of primes in Z Let us 


represent them as follows : 
news 

Ss, oe Kee teetrcoan 

Now let us take a look at the matrix 


diag (1,.--+¢1sS,41e---22 1) 


with s, on the position (k,k) and 1's in the other diagonal 


k 
positions. It is obviously that such a matrix can be 


decomposed into products of prime matrices as follows : 
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diag (17-.-17S, 4 1s--071) 


= diag (VseeelrP eleeee) appears te times 
1 

® diag (1s-+-1sPy 47 17++41) appears Cy times 
2 D 

® diag(1,--,1¢Pp «lee 1) appears r times. 
Pic nk 


However, equality (2) implies that 
A= U © diag(S,-1s-+-,1) 
@ diag (1,5, 71,-++01) 


e diag (1,---,1/S_) ® Vie 


which together with (3) yields the decomposition 


SiR ag Aas a (4) 
where 

PL = diag(1s--rlsPyeleeeeel) 
and Pie Pad are not necessarily different. As mentioned 
before, we can put in any unimodular matrices W, 
W-t between P. and Pe to obtain another decomposition. 


Given the decomposition (4), however, we choose instead to 


consider only 
Am (ire bee Ur mys (iene Us) ©... 9) (0 « ae e u-t) 
e (Ue Vv). 


Lemma 3.2. The matrices U ¢ P,e U-t are commuting normal 
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prime matrices. 
Proof : 
For each Pie there is a unimodular matrix 3 ' such 
that 
Gated P Vettel st = diagel, oe teh.) 


Let 


U, © (U-1 «© P, e U) © (U,)-1 
= (U,.' et « U-1) @ ae e (U « J-1 e (Up is") 
aoe e P.. e i ie 
= diag (1,.--.-1,p,)- 

Thus, each JU e Pie U-t ais a normal prime matrix. 


The commutativity is obvious. 


Q.E.D. 


We have therefore proved 
Theorem 3.3. Any nonsingular matrix A, which is not 
unimodular, in M(n,Z) has a left normal decomposition 
Dia 20) * oP er Uae £4. CU 5e Pe ere dee neue be OS) 
SOUT OU) ve 
where each P. is of the form diag(1,---71,p,515-++41) and 


k 


Py is a prime element in Z. 


3.4 : Unigueness of Decomposition 


For any nonsingular matrix A in M(n,Z), by some 
algorithm, we can transform it into its Smith normal forn 


A=Ue diag (S},7----+-+,S,) eV, 
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and then obtain a left normal decomposition 
A= (U © P,e U-!) © (U © Poe Ul) o.... © (U © Po Un?) 


°e (Uevy. 


Although the Smith normal form of a matrix is unique, 
there are many different pairs of unimodular matrices U, V 
such that Ue Ae V is the Smith normal form of A (e.g., 


different algorithms may result in different U, V). 


Example 
Let 
r 14 24 4 
ABe 4 | 
t 10 184. 
By taking 
Cem Least 4 
U=| | 
& 4 ae a 
© 2 -2 4 
v= | | 
L -1 ya ie 
we get 
YRepate ¥ 
em 2) 20%, 
= | | 
CAAOE (6) )4 2 


On the other hand, if 
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and 


u 0 6 J 


also results in the Smith normal forn. 


There is, however, a relationship between different 
pairs of matrices U, V. 
Lemma 3.4. Let A be a nonsingular matrix, which is not 
unimodular, in M(n,Z). 
rr 
A = (U.¢ P,* U-!) © (U © Pie U-1) e.... © (U © PL © 0-1) 
SUF OaN) 
and 
AS (UN oP e CU) St) ce (Gs SF Poe: (UN es 
a (01 «Doe (UNH t) ee (US ey 


be two left normal decompositions of A. Then, there exists a 


unimodular matrix X such that 


Uf e Vi = €l/detGs)) «Xe (0 *¢ Vp ¢ Ate X~ fe 4, 
where the matrix X satisfies the condition 

At 0 X @. A =.0 mod ({det (A) |) 
with congruence being elementwise congruence. 


Proof 3: 
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We have that 
(U' © U-1) oe A 
= Ut e { u-1 2e(U e Pie U-1) j e...9(U «¢ pr et-1) 
@ (U ro) V) 


Ut e (p,@e--*D. 2 Vv‘) e (Vtjy-1 e@ U-1 e (U e V) 


U' e [Pp }e---9P, os (U8) stFe 1a tse VX] 


e [(v')-t © (U')-1 © Ut] e U-t © (Ue Vy) 


Gee tee Se eye Shera. OU ei eae (0) —* | 
e (U' © Vt) © (Ut © Vty-t o Of e Q-1 
= (Ue) 
= Ae (Ut @ Vty-1 © (U' © U-t) © (U eV). 
This implies that 
(U' @ U-1) © A 
Sp Cer CUS Se ier See one Ue tee) Cle Vy s  (5) 
By taking X = Ut e U-! (hence X is unimodular) and by 
multiplying both sides of equality (5) by the 


adjoint, At, of A, it follows that 


At e2@X e@aA 

= det(A) * (U' e V')-1 e X e (Ue V). (6) 
Therefore, 

At  e°X © Arss0 mod (j{det (A) |) 
and 


ut = X e U. 
From equality (6) we have 

1/det(A) © At * Xe A 

o> (Um Ny, SS Oke ache Yr (7) 
But it can easily be shown that 


(V/det.cay) en At = x «kya? 
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1/det (A) e A+ @ X-1 eo A 
= 0 mod (j{det (A) |). 
Taking the inverse on both sides of equality (7), we 
finally obtain 
(1/det (A) © At © X e@ A)-1 
= ({ (U' e V'y-1 @ X © (U © V)}-1; 
that is, 
1/det (A) © At e X¥-1t oe A 


= (U0 e V)-! © X-1 e (Ut e Vt). 


U' oe Vt = 1/det(A) © (U ® V) © At e X-1 © A. 


Q.E.D. 


Definition 3.4. For any nonsingular matrix A in 
M(n,Z), a unimodular matrix X in M(n,Z) such that 
A+ e@ X @ A =O mod (| det (A) {) 


is called a unicongruential matrix of A. 


Proposition 3.5. For any nonsingular matrix A, which is 
not unimodular, in. M(n,Z), the set of all unicongruential 
matrices in M(n,Z) of A forms a subgroup, the 
unicongruential subgroup of matrix A , of the group of units 
-in M(n,2Z) {i-e., the set of all unimodular matrices in 


M(n,Z))- 


Let u be any nonzero element in 2, then the principal 


congruence subgroup of the group of units of level u is the 
set of all unimodular matrices W such that 
Wet mod (u). 


Further studies on this kind of subgroup can be found in 
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{Newman; 1972, chapter VII}. 


Proposition 3.6. Let A be a nonsingular matrix, which 
is not unimodular, in M(n,Z). Then the principal congruence 
subgroup of level det(A) is a subgroup of the uncongruential 


subgroup of the matrix A. 


From the propositions above we immediately get that the 
unicongruential subgroup of any nonsingular and non- 
unimodular matrix in M(n,Z) is a nontrivial group (a trivial 


group consists of the identity matrix only). 


Theorem 3.7. (Uniqueness of Left Normal Decomposition). 


Let A be a nonSingular matrix, which is not unimodular, in 
M(n,Z). Then, there exists a left normal decomposition 

A = Pact Eo eve ® Bae U. 
Moreover, if 

A = Q)° Qe ...9 O° V 
is another left normal decomposition. Then © = s, and there 
exists a matrix X in the unicongruential subgroup of A such 


that 
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Let A be a nonsingular matrix, which is not 
unimodular, in M(n,Z)- The existence of a left normal 
decomposition was proved in the section 3.3. Now 
assume that matrix A has another left normal 


decomposition given by 
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n= eG Merman Se Okey. 
Then by the fact that the determinants of matrices 
P. and Q, are prime numbers, it is easy to see that 
the number, s, of matrices Py should be equal to the 
number, vr, of matrices Q, es 

Secondly, by the fact that the matrices P, are 
commuting normal prime matrices, there exists a 
uNntmodulary) Batcix’ W so that Wot « P. © W = P*,, 
all i, where Pt is a diagonal matrix (see Thm 4.3 
for the proof of this fact). Similarly, there exists 
a unimodular matrix Y so that Y-1 e Q; e Y= Q's for 
all j, where QO", is a diagonal matrix. Thus, the 
decompositions become 
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We U' =U 
and 

Ye Vt = °¥. 


Furthermore, by a suitable rearrangement of the 
factors in the equations (8) and (9) and by Lemma 
3.4, we have 
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where X iS a unicongruential matrix of A. In 
addition, 


VY = Yvrex«' 


(1/fdet(A)) © X e (We Ut) © (At @ X-2 e A) 
=X @ Te (A-t @1X-1 © Ap. 


O.£.D. 


In concluding this chapter, we point out that the one- 
sided decomposition of matrices given above does not bear 
exactly the same meaning as_ the theorem of unique 
factorization in principle ideal domains, because we have 
restricted our 'prime factors' to be a subclass of all prime 


elements in M(n,Z). 
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CHAPTER IV 


Normal Prime Matrices in M(n,Z) 


The normal prime matrices introduced in the previous 
chapter possess some special properties, and we devote this 
chapter to exploring these. In section 4.2 we show that two 
commutative normal prime matrices can be diagonalized 
Simultaneously. The first section gives two lemmas which 


permit us to prove this result. 


4.1 : Preliminary Lemmas 


Lemma 4.1. Let P diag (1,---¢1¢P)» 


Q = U-1 © diag(1,---1,q) © U 
be two normal prime matrices with 
r Un-1 p14, 
v=" | l 
t v2 ee ee 
If Pe QgQ=QqepP andu #0 


then we have 
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with Q1 : (n-1)x(n-1) matrix 
Qt : (n-1) -column vector 
Q2 s (n-1)-row vector. 
Then, by assumption 
Pee vQe= One: Be 
so that 
e7Qh-t "Of, Ole. Oe pe 
1 =e ! 
t pQ2 px 4 tl Q2 xp J. 
But p > 1, so we get 
r 04 


and 

Q2 = (0,....-,0). 
Furthermore, from 

Q = U-1 e diag(1,-.-.-1,qg) ° U 
we have 


U e Q = diag (1s... 24/174) : U, 
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Then we get the following equations ; 


UuUn-1 e Qn-t = Un-1, 
Vie xe = Ul, 
ux = que 


But by assumption u # 0. Thus from equation 
follows that 

x = q # 0. 
Combining this with equation (2) we have 


Ut @ x = [Jt e g = U1, 


so that 
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dae | 
to Qian 


Moreover, because UJ is unimodular then u = +1. 
Therefore, from 
u ¢ (det (U"—*)) 
= +(det (Un-!)) 
= det (3) 
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We have proved that U"-1 is a unimodular matrix. 
Now by equation (1), we conclude that 

Qn-1 =} 
which completes the proof. 


Q.E.D. 


Lemma 4.2. Let P diag (V7.2 71 DP) ¢ 
OF= 0-1 (evdiag ll ca. ¢ 1 ag) oe U 
be two normal prime matrices with 
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with Q2-1 being a normal prime matrix. 

Proof : 
As in the proof of Lemma 4.1, the matrix Q must be of 
the form 


POU. On 


Furthermore, U! # O, otherwise U is not unimodular. 


Now from 


Ue Q = diag(1,---.41,9) ° U, 
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and we get the following equations : 
UM-1 e Qn-1 = Yn-1, (4) 
U1 ex = 1. (53 
From equation (5) and the fact that U! # O we get 
xX -Spt; 
and therefore 
[det (O2-*y | = (det (0)! = q- 
Now, in order to prove Q"-1 is normal prime we need 
to find a unimodular matrix T"-1 such that 
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But, if we let 
Pena tho 4 


Te= 4 | (6) 


twee (mies We 
Le Pan de 


Pe | od ae 
' whl 7 x3 . 
at hs fey ze “¢ ay ty pre y. 7 ; 
j sae " ; j P 
1 © - 7 al al, 7 
a ' Fit i ‘. 7 : D 
us oot! Son See ee ' 
~ { oi i A 17 
fi if £ “TE iy 
' : ; , : a: 
: re ~ 
! / a) P 
ft y.07 2 b. 
. 5 ' } } | a % t+! i 
i > | } 
» | : ‘ } iS, 'e SO 
‘9 pnbwotiot sit teh Se hae 
gi-@y a. $83.05 29 ) 
5 in‘ ott 2 g e yy ~ :* 
é : | a 
i] 7 
w OR 1 Cee One See ee: ey, notshupy wots 
: ~ ia — bp & 
_ a 1 J 5 7 ; : v4 ; 
» ‘a= Pie Ans aF 7's ath 
4 - ia . ey =o 
rat fy oll oa 7 
yt aes seria han : 
P me 6 
' E ras 
v= ay : 
; aR ‘= esa ee Jattooptebi f 


at a amy Te | vaduai rp ae “Sapp, Pee ce enae i 8 
it Cond ye: a2 i" 


Aine on aaa ene a Bet? os 
: | Wa 
Tey, o. une — 7 is 


"i ir 7 


a Z ey FY 
ee ee 


33 


Ce ee Ober OC eee eOeg es (TO-t)—t 0 4 


i) 


diag (1, «<¢ 1s9,¢ 1) 


Y-t e diag(1,---41,q) © Y, 


with 
Re TPR OO <4 
Y= |] =O Oe & 
ve) 8) Th (Oras 

ives 


GY eet hy esG. sO (Yregt)—* 
= diag (1, <i, 174) < 
Thus, if there exists a matrix T"-1 which transforms 
Q®-1 to its Smith normal form, then the matrix V = Y 
e 7 transforms Q into its Smith normal form. The 
matrix V = Y e T assumes the Special form 
p Via “vi. O54 


Yas {°°O 0 Vel (7) 


Therefore if such a matrix V can be constructed, then 
T2—-1 can be obtained from T = Y-! e V and will assume 
the form (6). 
From equation (4) and from the fact that [det(Q —!) | 
= q, we get that U"-! must be a singular matrix. 
Moreover, by Theorem 1.2, there exist two unimodular 
matrices H[-1 and K"™-! such that 

HoO—-1 »® U-1 e KN-!1 


= diag(ut,... ,u 2,0) 
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Therefore, let V=S e (H e¢ U), and this is the 
desired unimodular matrix of the form (7). 


Q.E.D. 


4.2 : Diagonalization of Normal Prime Matrices 


SS eee 


Any two matrices A, B in M(n,Z), even if they are 
commutative, can not in general be diagonalized by the same 
sequence of elementary operations, The situation, however, 
is quite different if A and B are both normal prime 
matrices. 


Theorem 4,3. Let P, Q be two normal prime matrices in 
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Pree = =O) 6 P 
then there exists a unimodular matrix U in M(n,Z) such 


that U-1 « Pe UY and U-t« Qoe 9 are both of diagonal 


Step 1: 
By the definition of a normal prime matrix, 
there is a unimodular matrix U' such that 
Oa tresP ser UO = "diag (Ven cay Ve Bs 
Step? 203 
Apply Ut to the matrix Q to obtain 
U'-1 © Qev= Qt. 
Since: 40% is still commutative with 
diag(1,--.-,1,p), by lemma 4.1 and lemma 4.2, Q! 
can assume one of two forms. 
Step 3% 


rf 


ele: qf, 

then Ut is as desired and we have finished. 
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(UNn—1)-4 e g'n—-i1 e gun-1 


= diag(1,.-.-.,1,q).- 


Step 5: 
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By induction we can easily generalize the theorem above 
for m commutative normal prime matrices. 


Theorem 4.4. Let Pt, P2,...,P" be normal prime matrices. 


If they are commutative, then there exists a unimodular 


matrix U such that U-1!1 e pte 9 is a diagonal matrix for 
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In concluding this chapter we give two examples to 
illustrate Theorem 4.3. 
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Hence, P and Q are commutative. 
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The matrix U is then given by 
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CHAPTER V 


A Non-iterative Algorithm for Computing 


A Left Greatest Common Divisor 


5.1 : The Algorithm 


Unlike the case of integers and polynomials, because of 
the non-commutativity of matrix rings, the prime 
decomposition of two matrices does not yield their left 
greatest common divisor. Having defined an algorithm for the 
division of two matrices, however, the Euclidean algorithm 
can be used to obtain a greatest common divisor. In this 
chapter, a new algorithm for the computation of a left 
greatest common divisor is given. This algorithm requires 
only the computation of the Smith normal form of a matrix, 
and thereby avoids the iterative nature of the Euclidean 


algorithm. 


We begin with a definition of a left greatest common 
divisor, more general than that given in Definition 3.2. 


Definition 5.1 : Let A, B be arbitrary elements in M(n,z), A 


left greatest common divisor of A and B is an element D, 
denoted by lgcd(A,B) = D, in M(n,z) such that D left divides 
both A and B; and if D' is any element in M(n,z) which left 
divides both A and B also, then D! left divides D. 


The Non-iterative LGCD Algorithm : 


Given two matrices A and B in M(n,z), this new 
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algorithm computes D = lgcd(A,B) and two multipliers X and Y 
such that 
Be 3 SFB ve Ky 7 Ds 


Step 1 (Smith normal form) : 


i] 


Augment the matrices A, B to form an nx2n matrix C 
{ A,B Jj. Then find two unimodular matrices U, V 
such that U eC *eV= [5S ,0 ]is the Smith normal 
fornsof C. 

Step 2 (Inversion) : 
Compute U-?. 


Step 3 (Termination) : 


D <--- U-1 e (the first n columns of U @ C @ Vj). 
X <--- the first n rows and first n columns of V. 
Y <--- the last n rows and first n columns of VY. 


roof of the validity of the Algorithm : 
By Theorem 1.2, there are two unimodular matrices 
eneayye V(2nx2n) such that 

US Cle yess iSa,-0: ] 


is the Smith normal form of C. Partition the matrix V 


into four nxn submatrices as follows :; 


Hence, 
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Ave V@.+ Be Ys =D. 
To prove D is a left greatest common divisor we show 
that D left divides both A, B and any other left 
common divisor of A, B left divides D. The second 
part is quite obvious, and we need only to prove the 
first past. 

Since the matrix V is unimodular, its inverse 
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where Wi are nxn matrices, exists. 
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A=DeW, 
B= De We. 
Thus, D left divides both A and B. Furthermore, by 
setting X = V1, Y = V3, it follows that 


A eX + Be ¥Y = PD. 


Q.E.D. 
5.2 : The Algorithm for Multiple Matrices 
The previous algorithms can be applied m - 1 times to 


compute a left greatest commom divisor of m (m > 2) matrices 
in M(n,Z). It turns out that by making a small modification 
to the non-iterative LGCD algorithm we get a new efficient 
algorithm for computing their left greatest commom divisor. 


The Non-iterative Algorithm for m Matrices : 


Given: m, (@ >e2) “matrices Al, ‘A220, 5A" in 4 (0;2).> 
this algorithm computes D = lgcd (At,A2,..-,A®%) and m 
multinliérs: X%*, %2, 20, X" such that 
AAS 10 XA ASE itil eo X= D 
Step 1 (Smith normal form) : 
Augment the matrices At, A2, ..., A™ to form an 
nx (mxn) matrix C = [{ At , A2 , ..- , A™ jj Then find 
two unimodular matrices U, V such that 
DreeCiesVesesius.p ON, Noes F OM] Is: thevvSaLlthytnornal 
form of C. 
Step 2 (Inversion) : 
Compute U-!. 


Step 3 (Termination) : 
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D <--- U-1 e (the first n columns of U e C ® V). 
X1 <--- the first n rows and first n columns of VY. 
X2 <--- the second n rows and first n columns of V. 
x™ <--- the last n rows and first n columns of V. 


5-3 : Complexity Considerations 


There are now three methods of computing a greatest 
common divisor lgcd(A,B) of two matrices A and B, namely, 
Method (1) : 
Euclidean algorithm using Sanov's division algorithm 
(see Section 2.2). 
Method (2) : 
Euclidean algorithm using the integer-arithmetic 
algorithm (see Section 2.3). 
and 
Method (3) : 
The non-iterative lgcd algorithm given in the 


previous section. 


We first count the number of operations required to 
perform the divisions in methods 1 and 2. Sanov's algorithn, 
in steps 1, 2 and 3 includes the triangularization of a 
matrix, a matrix inversion and a matrix multiplication. This 
requires O(n3) operations. Steps 5 and 6 require an 


additional 0(n3) operations. Thus, Sanov's division 
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algorithm requires at least O(n3) operations. 


On the other hand, the integer-arithmetic division 
algorithm involves a matrix multiplication, the computation 
of the adjoint of a matrix, and finding the Smith normal 
form of a matrix. The Smith normal form of a matrix can be 
obtained by means of Bradley's algorithm, which requires 
O(n3) + O(n2 © determinant of the matrix) operations 
{Bradley; 1971]. Indeed, this is the dominating cost of the 
integer-arithmetic division algorithm. Thus, with respect to 
total operations required, this crude analysis does not 
permit us to determine which of the methods 1 or 2 is 


superior. 


Methods 1 and 2 are iterative, requiring at most 
min{ j|det(A){, {det(B){| } steps (i.-e., matrix divisions) 
before termination occurs. Method 3, on other hand, requires 
Simply one matrix inversion, and the Smith normal form of 
one nx2n matrix. Clearly then, method 3 is asymptotically 
superior to both methods 1 and 2 with respect to the total 


number of operations required. 


The above analysis has ignored the cost of each of the 
operations involved. Method 1, for example, uses rational 
arithmetic; and even if matrices A, B have single-precision 
components, all methods may require multiple-precision 
arithmetic. Thus, the cost of the methods depends not only 
on the total number of operations required, but also on the 
size of intermediate results on which multiple-precision 


arithmetic is being performed. Indeed, we suspect that all 
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three methods suffer because of “intermediate expression 
growth", a phenomenon common to most algorithms in algebraic 
and symbolic manipulation. That is, the length of 
intermediate results ae be large even for problems where 
the length of the initial and final results is small. An 
analysis of the three methods which takes these matters into 
consideration, however, is a major undertaking, and we leave 


it as a subject for further research. 
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